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' We implement in systems of fermions the formalism of pseudoclassical paths 

o ■ 

CN ■ that we recently developed for systems of bosons and show that quantum states of 

Qh. fermionic fields can be described, in the Heisenberg picture, as linear combinations of 
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randomly distributed paths that do not interfere between themselves and obey clas- 
sical Dirac equations. Every physical observable is assigned a time-dependent value 
on each path in a way that respects the anticommutative algebra between quantum 
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operators and we observe that these values on paths do not necessarily satisfy the 
O ' usual algebraic relations between classical observables. We use these pseudoclassical 

■ paths to define the dynamics of quantum fluctuations in systems of fermions and 

show that, as we found for systems of bosons, the dynamics of fluctuations of a wide 
Qj I class of observables that we call "collective" observables can be approximately de- 



scribed in terms of classical stochastic concepts. Finally, we apply this formalism to 



^ . describe the dynamics of local fluctuations of globally conserved fermion numbers. 

PACS numbers: 11.10-z, 02.50Ey, 03.65Yz 



I. INTRODUCTION 

Quantum fluctuations of physical observables can be manifest when they couple to an- 
other external system. For example, the Lamb shift in the absortion/emission lines of atoms 
is a direct consequence of their interaction with vacuum fluctuations of the electromagnetic 
fields. Similarly, quantum fluctuations in the energy-momentum tensor can affect the space- 
time curvature and could result in dissipative effects in quantum cosmology and quantum 
gravity 

The fluctuations and their dynamics can be described using the general formalism of 
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time dependent Green functions to obtain the correlation of physical observables at different 
times. Some aspects of quantum fluctuations of commuting physical observables, like photon 
shot noise in optical systems or the Casimir force between plates, can be simulated as classical 
stochastic fluctuations with well defined probability distribution functions. However, it is 
not clear in this context to what extent it is possible to describe additional aspects of the 
dynamics of quantum fluctuations using classical stochastic concepts. 

In bosonic systems classicality is commonly associated with highly populated coherent 
states that can result, for example, from high temperatures or strong particle production 
in cosmology. In this context issues concerning classicality and decoherence of quantum 

nn nn 

fluctuations p|, \^ during and after inflation, and their time evolution p, la| were addressed. 
Fermionic statistics, on the contrary, prevents any state to be highly populated and, more- 

n 

over, fermionic coherent states cannot be identified with classical trajectories [21, so it is 
generally conceded that there is no issue of classicality in quantum fluctuations of fermionic 
fields. 

In a recent paper we have developed a new formalism of pseudoclassical incoherent 

paths to describe the dynamics of quantum fluctuations. New aspects of classical behaviour, 

in particular classical dynamics and decoherence, emerge in this formalism associated to 

fluctuations of collective observables which depend on a large number of degrees of freedom, 

even in scarcely populated states like the vacuum (a conclusion anticipated in 0] by different 

means). The notion of classicality that we propose is quite different than the case which is 

often discussed in the literature: it is associated to collective observables rather than to a 

n 

certain class of quantum states, namely coherent states. In ^8] we elaborated these ideas for a 
system of weakly interacting bosons. In this paper we show that this notion of classicality is 
not exclusive to systems of bosons, as it also emerges associated to fluctuations of collective 
observables in systems with a large number of fermions. 

In the new formalism quantum states can be represented, in the Heisenberg picture, 
as linear combinations of randomly distributed pseudoclassical paths (PGP's) that do not 
interfere between themselves. We use each set of random PGP's to deflne the dynamics 
of quantum fluctuations in that state of the system. We remark that any single PGP 
cannot be identifled with any quantum state. Instead, a whole set of PGP's with their 
corresponding probabilities to randomly happen at any time is the only formal object that 
can describe the dynamics of quantum fluctuations in one state of the system. On each 
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path of a set every physical observable is assigned a time-dependent value in a way that is 
consistent with the commutation/anticommutation relations between quantum operators. In 
particular, the generalized canonical coordinates and their conjugate momentum operators 
get time-dependent values which obey classical equations of motion and depict a collection of 
harmonic oscillators with constrained initial conditions. On the other hand, we notice that 
these pseudoclassical time-dependent values on paths do not necessarily respect the usual 
algebraic relations that classical observables fulfill, as a consequence of the non-commutative 
relations between quantum operators. 

We select collective observables which depend on a large number of bosonic or fermionic 
independent degrees of freedom because we realize that their time-dependent values on 
paths do approximately regain the usual algebraic relations between classical observables 
and, therefore, we can use the formalism of PGP's to describe the dynamics of quantum 
fluctuations of collective observables in terms of unconstrained classical stochastic processes. 

The formalism of PGP's, for bosons or fermions, shows that we can give a description of 
quantum mechanical states in terms of non-interfering paths if we trade it for a non-trivial 
definition of the algebraic relations between the values on paths of physical observables. This 
formalism can help to understand the process of decoherence and the onset of classicality in 
quantum systems because it can be applied to closed systems without explicit reference to 
observers, measurement, or an environment. In some cases it may also become a useful tool 
for performing calculations, in particular numerical simulations. 

This paper is organized in five sections. In section II we review the fundamental concepts 
of the formalism of PGP's and implement them in a system with a single fermion mode. In 
section 111 we extend the formalism to a system with two non-interacting fermions, and then 
we generalize it to systems with any finite number of free fermions. In section IV we define 
collective observables and describe their dynamics in the context of a regularized QFT. 
Section V contains a summary of results. This paper deals basically with linear fermionic 
fields. A discussion of the general case of interacting fields is postponed to a forthcoming 
work. 
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II. PSEUDOCLASSICLAL PATHS OF A SINGLE FERMIONIC MODE 

In this section we present a detailed description of a system with a single fermionic mode 
in the new formalism of PCP's. We will start building the set of random paths which 
describes the dynamics of quantum fluctuations in the vacuum state and, afterwards, will 
show how to build the sets of paths that describe states other than the vacuum. 

The hamiltonian of the system is 

H = n (a^a - I) , (1) 



2. 

where the operators a"^ , a, which create and annihilate the excitations of the mode, obey 
anticommutation relations {a^,a} = 1 and {a"'',a^} = {a, a} = 0. The last two relations 
imply, in particular, that (a^)^ = = and, therefore, the orthonormal basis of eigenstates 
of the hamiltonian, which linearly span the Hilbert space of the states of the system, contains 
only two independent vectors: {|0), |1)}. The action of the operators a\ a on this basis is 
given by the equations a\0) = 0, a\l) = |0) and a^O) = = 0. We will fix for 

simplicity k = 1 and, therefore, H\0) = — 1|0) and H\l) = 

The operator is the hermitic conjugate of the operator a, so that the couple of operators 

e = ^(« + «^), e = -^(a-at), (2) 

are hermitic and obey anticommutation relations: {■CjO = {'CjO = 1 ^"^^ — 0- 

The action of these operators in the basis of eigenstates of the hamiltonian is summarized 

in the equations ^|0) = = ^|0) and i\0) = = ^\0), and, therefore, 

they can be identified with Pauli matrices, ^ = :^cri and ^ = ~:^cr2- Notice that = = 
1 

2" 

The most general linear operator that can be defined on the Hilbert space of the single 
fermionic mode is a polynom of the type, O = ai ■ 1 + 02 ■ ^ + «3 ■ ^ + «4 ■ i^^, where the 
coefficients aj,j = 1,2,3,4 are, in general, complex numbers. This linear set of operators 
is the Grassmann algebra generated by the anticommuting hermitic operators The 
operator C, (or independently, ^) forms a complete representation of commuting observables 
in the Hilbert space. An operator O is hermitic if and only if the coefficients aj,j = 1, 2, 3, 4 
are all real. In particular, the hamiltonian is 

H = ii^ (3) 
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Our first aim is now to write the action of a generic operator O on tlie vacuum 
state in terms of the action on the same state of some linear combination of the iden- 
tity operator 1 and the generator ^. In general, it can be immediately checked that 
O|0) = {^{ai — \a4) ■ 1 + {a2 — ia^)^^ |0) = Vo{l, is also immediate to check that 

the operator 7^0(1, defined through this equation is unique as the identity /3i-l|0) = /92-^|0) 
only holds if (3i = (32 = 0. From the identities, 

m = -^m, H\0) = -l\0), (4) 

we obtain, for example, that V^l,^) = —iC, and Vh{l,0 = ~i ' ^■ 

The next stage in our programme is to identify the two eigenstates of the operator ^: 
= -^(|0) + |1)) and = -^(|0) — |1)), and their corresponding eigenvalues, ^+ = ^ 
and ^- = Then, we expand the vacuum state in the new basis: |0) = + 
= J2q=qj^4'{(l)\(l) , with 4'{q±) = The wavefunction of the vacuum in the basis of 
eigenstates of ^, ip{q), is equally valued on each of the two disconnected points \q+), \q-) 
of the configuration space and so the probability \4'{q±)\'^ = (:^)^ = I each of them to 
happen randomly. 

Moreover, using the relation ^|0) = ^(^|g+) + ^1^-)) = :^(A+|g+) + A_|g_)) we can 
assign to the operator ^ the classical value 

on the first point, and 

on the second point of the configuration space. Thus, we have defined on configuration space 
a random variable C,ci{q)- 

Furthermore, following the identity O|0) = Po(l,0|0) = 75(^0(1, + 
Toil, C,ci{q~))\q-)) , the generic operator O should be assigned the random variable 

Oci{q±) = Vo{l,U{q±)) = -, — r^T , (7) 

that gives the operator a pseudoclassical value at each point of the vacuum configuration 
space. For example, looking at Q we assign to the operator ^ the random variable C,ci{q) 
which takes the value 

— — z 

^d{q+) = ~i^d{q+) = ^ (8) 
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on the first point of tlie configuration space and 

liiq-) = -^Uq~) = (9) 

on the second point. Also according to the free hamiltonian Q is assigned on each of 
these points the constant value 

hM±) = -l = m\o). (10) 

Let then go a step further to describe the time dependence of the random variables 
we have just defined. In the Heisenberg picture, time dependence of the operator O is 
described by the expression 0(t) = e*^*Oe~*^*, where H is the hamiltonian of the system. 
The observable 0{t) can then be associated following the same algorithm that we have 
described above with a new random variable Oci{t,q±). Thus, each physical observable is 
given a time- dependent c- value at each point \q) of the configuration space, that is a path 
whose realization probability is |'?/'(g)p. 

In particular, the time evolution of the operators ^ and ^ is described by the expressions 

m = e'^^'ic-'"' = cos{t)i - sin{t)l (11) 

lit) = e'^'ie-'^' = sin{t)^ + cos{t)i, (12) 
which solve the differential equations 

^-i,, ^-m. (X3) 

On the vacuum state |0) these equations imply (see (jl}) 

e(t)|0) = e^*e|0), e(t)|0) = -^e^*e|0), (14) 

which, according to our formalism, mean that in the ground state |0) the time-dependent 
classical values of this pair of operators on each of the two paths |g-|-) are ^c/(^, <i±) = e**^cz('?±) 
and iciii^ 9±) = ici{.(l±) ■ They obey the classical equations of an harmonic oscillation 

= -icAt] g±), = icAt] g±), (15) 

with initial conditions fixed by (0 and |HI) on the first path, and (0 and E)) on the second 
path. The operator i^^ commutes with the hamiltonian so it does not change with time 
and, therefore, neither its classical value (fTUIl on each path does. 
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The formalism of PCP's can be extended to describe any other state in the Hilbert 
space of the single fermionic mode. Consider, for example, the normalized state |\E') = 
cos{9)\0) + sm(6')e*'^|l) and let expand it in the basis {|g'±)} of eigenstates of the generator 

1^) = -L ((cos(^) + sin{9)e'''')\q+) + {cos{9) - sin{9)e'^)\q^)) . (16) 

On each of the two points in the configuration space the random variable C,ci{(l±) 
takes the value ^^^^j^^ = identical to those specified by eq. (0) and (P)- In order to 
obtain in the new quantum state the random variable which corresponds to the generic 
operator O we must find in the linear subspace spanned by the identity operator 1 and the 
generator ^ the operator ^^(1,^) such that 

O\^) = Vf{l,0\^)- (17) 

This can be done through the following three steps: 1) we write the state 
as a linear combination of the operators 1 and C, acting on the vacuum: |\E') = 
[cos{9) ■ 1 + V2sin{9)e''t'^) |0); 2) now, use the rules that we have stated before in this 
section in order to obtain the action of the operator O on the state |\E') in terms of the 
action of some new linear combination of the operators 1 and ^ on the vacuum state: 
0|^) = O (^cos{9) + V2sin{9)e'^^) |0) = ^(1,0 10); 3) note that the relation introduced in 
step 1) can be inverted into the relation |0) = cos^e)-sin^{e)e'^''t- {'^os{9) - y/2sin{9)e"^^^ |\E'), 
which is then introduced in the expression that we obtained from step 2) in order to get, as 
desired, the action of operator O on state |\E') in terms of the action of some linear combi- 
nation Vf{l,^) of the operators 1 and ^ on such quantum state. The three described steps 
of this algorithm are well and uniquely defined, so the operator Vf{l,^) is unique for each 
operator O. 

In particular, for the operator ^ we obtain the identity 

m = ^^^2^0^ _]^^2^0y., i^MO)cos{9)e^^ - {cos\9) + .zn^(^)e^^^)OI^), (18) 



which means that this operator is assigned, when the system is in the quantum state |^), 

(g±l*> 



the random variable i{q±) = ^l^^f^) whose value at each point of the configuration space is: 
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The dynamical equations (fTSjl . when apphed on the quantum state 

(20) 

enforce that classical values of the operators ^ and ^ on each pseudoclassical path of the 
random set that describes the state |\E') still evolve according to classical equations of motion 
(fT^. Initial conditions are now fixed on each path by ^ and ©, respectively, for C,ci{t = 0) 
and (|T9|l for C,ci(t = 0), instead of (jH)), In general, the state of the system fixes only 
the initial conditions for its set of PCP's, while the dynamics of these paths subsequently 
follows classical equations. 

The action of the hamiltonian (jS)) on the quantum state |\E') can be expressed, using (fTHjl . 

as 

(21) 

which implies that the two paths of the set that describes quantum fluctuations in the generic 
state |\1') do not necessarily have the same energy, but 

The two pseudoclassical paths either are no longer equally probable as can be seen from 
equation (fTHjl : when the system is in the quantum state the probability of the point 
to happen randomly is |?/'(g+)p = ^\cos{9) + sm(6')e*'^p, while the probability of the point 
is \ipiq-)\^ = ^\cos{e) - sin{e)e"^\'^. 

Notice that although expressions (fT^ and (f^^ are ill-defined at the event when 
g«(/) _ _|_]^ g^]^^ sin{9) = ±cos{6), the probability of such event to happen randomly is zero, 
while the other event \q+), then, should have an absolute probability equal to 1 to happen. 
If we consider this case as a limiting case when 9 — > we can see that event 

becomes less and less probable, but the corresponding values for the variables id{q-) and 
hci{q-) grow towards infinity. In the case of the latter variable the probability of the event 
decreases faster than the divergence of the variable, so that this event contributes and 
infinitesimal amount to the average or expected value of the classical value hd over the set 
of paths. On the contrary, the divergence in the expression for the former variable grows 
fast enough to contribute a finite amount of half the value of its total average < ^ > over 
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paths, even though the event has a tiny probabihty to happen. The situation is similar for 
the event when e*"^ = ±1 and sin{9) = ^cos{6). 

Let us add that although ^ci{t,q±) and ^ci{t,q±) obey classical equation of an harmonic 
oscillation, the energy of the PCP, as defined by equation (j221), is not necessarily equal to 

mcir+iii)')- 



We will finish this section with two results that were already noticed in 



but we remind 



them here for their relevance in the discusion of the next sections. The first result: 

{Mj\0\^) = = 0,Kg+)|^(g+)P + 0,Kg-)|^(g-)P =< >pcp, (23) 

means not only that the average of pseudoclassical values of the observable O over the whole 
set or ensemble of random PCP's that describes the state is equal to the average value 
of the operator O in this state, but it also means that PCP's do not interfere between 
themselves. 

The second result on PCP's that we want to bring to attention is related to the second 
momentum of the distribution of pseudoclassical values of the generic operator O over the 
ensemble of random paths: 

< {O^OU >pcp= {^\0^0\^) = (vl>|(P„^(l,0)^Po^(l,OI^) =< 0:P,i >pcp . (24) 

It means that the average on PCP's of pseudoclassical values for the operator 0^0 is equal to 
the second momentum of the distribution of Od over the paths. Although this equation holds 
in average for any operator O it does not necessarily hold on each one of the paths as can be 
readily seen, for instance, from equation (fT^ : while {i^i)ci = (l^)d = ^, we find after some 

1 K i-l, i- fC \*fC \ — (l+5S*«^(2f)(l+cos(2(/-))=F2sm(26»)(cos2(e)cos(0)+sm2(6»)cos(3</>))) 

algebra tliat [U) - 2(i-ism2(2e){i+cos(2^))) ' ^ 

sign in the numerator stands for the value of the random variable at and 1^-.), respec- 
tively. Only after averaging these two possible values with their corresponding probabilities 
we recover the identity ()24|). This is the reason why we call these paths pseudoclassical 
paths: although they obey classical equations of motion, the values on paths of physical 
observables which do not commute with the representation of commuting observables we 
have chosen (in this example, the operator C,) do not fulfill the usual algebraic relations 
between classical observables. Or stated in other words, {0''0)ci — {Oci)*{Oci) is not neces- 
sarily zero on each one of the paths, neither even on any of them. This is a consequence 
of the non-commutative relations between operators in the algebra of quantum observables. 
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We will return to discuss this aspect of the formalism at the end of next section when we 
will introduce collective observables and show that classical algebraic relations are naturally 
recovered for this kind of operators. 

III. PSEUDOCLASSICLAL PATHS IN A SYSTEM OF MANY FERMIONS 

In this section we discuss how to extend the formalism of pseudoclassical paths to systems 
that contain more than one fermion mode, for example a fermionic quantum field. The basic 
ideas are the same that we have already presented in the previous section, but there are 
now some subtleties related to the anticommutation relations between operators associated 
to different fermionic modes that we would like to notice. We will start, for the sake of 
simplicity, with a system with two free fermions whose hamiltonian is: 

H = Ka (ata -l)+^b (b^b - = if„ + H,. (25) 

The operators a^,a create and annihilate excitations of a first fermion mode, and 6^,6 create 
and annihilate excitations of a second fermion mode. They obey anticommutation relations 
{a^,a} = {b\b} = 1, and the anticommutator between any other pair of these operators 
is equal to zero. In particular, those relations imply that (a^)^ = = (fe"!")^ = = 
and, moreover, the operations that create or annihilate excitations of the first fermionic 
mode anticommute, instead of commuting, with the operations that create or annihilate 
excitations of the second fermionic mode. 

The basis of eigenstates of the hamiltonian can be obtained applying the creation 
operators and 6^ on the vacuum state |0; 0). The anticommutation relations between the 
operators guarantee that the Fock space they generate contains only antisymmetric vectors: 

|0;0) |1;0) = a"^|0;0) |0; 1) = 6"^|0; 0) = a^b^\0;0). (26) 

Their corresponding eigenvalues can be read of the equalities 

i/|0;0) = -^^|0;0) i/|l;0) = +^^|1;0) (27) 

iJ|0;l) = -^^|0;1) = +^^|l;l) . (28) 

The operators a^, 6^ are the hermitic conjugates of a and b, respectively, so we can define 
two couples of hermitic linear combinations of them: the hermitic operators introduced 
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in (j21), together with the couple 
These four operators anticommute each other 

{e, o = {e, c} = {e, a = {I c} = {I (} = {c, a = o, (30) 

and each one is proportional to the identity operator when it is multiplied by itself, = 
= = = |. They generate the Grassmann algebra of linear operators defined on the 

Hilbert space spanned by (j26|l . 

The two components Ha and Hh of the hamiltonian (j25p can be expressed in terms of the 

new operators as Ha = i^ai^iO Hf, = Ubi'^COj which commute each other [Ha, Hi,] = 0. 

Therefore, e*^* = e'^^'^^^"^^ = e'"-^e'"^\ and 

^{t) = e'^'ie-'"' = e'^-^^ee"'^"*, (31) 

while 

((t) = e'^\e-'"' = e'^^Xe-'^'K (32) 
Similar expressions are obtained for ^{t) and Q{t) and they imply 

Kam (33) 

K.at). (34) 

In order to build the set of PGP's which describes a generic quantum state |^) of this sys- 
tem we need to choose a complete representation of commuting observables. The operators 
^ and C do not commute each other and, therefore, are not eligible as such a representation. 
Instead we use the pair of hermitic operators ^ and i^C, that, according to do commute 
each other. 

So we expand the state = '}2q'^{l)\'i) the basis of common eigenstates \q) to the 
complete representation of commuting observables that we have chosen. We know from 
the discussion of the previous section that each of these eigenstates is actually promoted 
to one PGP, whose probability to happen at any time is |\I'(g)p. Each observable of our 
representation is naturally given as its classical value at each of these eigenstates its own 



(6-6t). 



(29) 



dm 

dt 



dm 

dt 



dm 

dt 



dm 

dt 
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eigenvalue on it and then, exploiting identities of the kind 0(t)|\&) = Vf{l, iC,C', any 
other physical observable O is also assigned a time-dependent value on each of them: 

Oci{t;q) = Vtit;l,Uq),mciiq))- (35) 

We can see from this last expression that whatever the quantum state |\E') is, equations 
and (jH^ enforce when they are applied on such state (see eq. ((201)) that the two couples 
of observables ^, ^ and (, ( will be assigned following this formalism time-dependent values 
that obey classical equations of two decoupled harmonic oscillators, 

d^cl{t,q) r , . dicl{t,q) , . . . 
= -KaU{t, q), — = KaU{t, q) (36) 

= -KbCci{t, q), = KbCci{t, g), (37) 

on each PCP, generically labeled as of the random set which describes that state 

As an example we will explicitly construct here the set of PCP's that describes the vacuum 

state |0;0) of this system. The generalization to other states is straightforward following 

the three-steps algorithm discussed in the previous section and we will only outline how to 

apply it to the two-fermions system. 

The first stage then is to find the way to express the action on the vacuum state |0; 0) of 

any operator O in the Grassmann algebra generated by ( and ( in terms of the action on 

the same state of some linear combination of the set of operators generated by the complete 

representation of commuting observables, ^ and iC,C'- 

O|0; 0) = (ai ■ 1 + as ■ e + «3 ■ (^^C) + «4 ■ 0) = ^o(l, ^^OlO; 0). (38) 

The operator Vo{l,^,i^() defined by this expression is unique because the action of the 
identity operator 1 on the vacuum does not excite any mode, the action of the operator ^ 
excites the mode |1; 0), the action of the operator (^^C) excites |1; 1) and the action of (i^^C) 
excites |0; 1). Therefore, the only linear combination of these four operators which can give 
the number zero when acting on the vacuum is the trivial combination. 

For example, we already know that ^|0; 0) = — i^lO; 0) because of the identity a|0; 0) = 0. 
A little more tricky is the way we describe the action of the operator ( on the vacuum. We 
first note that C = 2^(f'C) = -'^{iOl Therefore, C|0;0) = -2(|'C)I^|0; 0), and subtituting 
in this expression the previous result we get C|0;0) = 2{i^Q^\0]0) = 2^{i^()\0; 0) . Now, 
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using the identity 6|0;0) = 0, we get C|0;0) = -<|0;0) = -2i^{i^C)\0; 0) . Using the 
anticommutation relations between the generators of the Grassmann algebra we can write 
down similar expression for any other operator O. 

The next stage is to obtain the basis of common eigenvectors to the complete represen- 
tation of commuting observables we have chosen,^ and («^C)- 

= I m 0) + |1; 0)) + (|0; 1) + |1; 1))) (39) 

g_) = ^ ((|0; 0) + |1; 0)) - (|0; 1) + |1; 1))) (40) 
|g-;g+) = ^((|0;0)-|l;0))-(|0;l)-|l;l))) (41) 

|g_;g_)^i((|0;0)-|l;0)) + (|0;l)-|l;l))), (42) 

and their corresponding eigenvalues, 

1 ~ 1 

^1?+; Q+) = 1+) ^^C\q+; q+) = +2 1^+- ('^^) 

Q-) = iickv, Q-) = -^1^+' (44) 

= iiC\q-; q+) = +^\q-; q+) (45) 

In the new basis fjS^ - iP^ the vacuum state can be expanded as 

|0;0) = i(|g+;g+) + |g+;g-) + |g_;g+) + |g_;g-))= ^ ^{ql,qi)\ql;qi). (47) 

The configuration space of this system has four disconnected points, each one of them has, 
when the system is in its vacuum, a probability |'?^(5')P = (|) to happen randomly. 

The classical values at each point of the configuration space of the operators that form our 
complete representation of commuting observables can be obtained directly from (j43|) - (j^6|) . 

For example, Ul^^q^) = ^^^^^ = ^ and mUqv,q^) = ^'^^^T = ^ ^'^^^^'^Y^ 
we define their classical values in the other three points of the configuration space. 

Now we can benefit from identity (jHHj) in order to assign to any other opera- 
tor O its classical value at each point of the configuration space: Oci{q±,q±) = 
'Po{l,^ci{q±,q±),{'iiOd{q±,q±)) = ^'(^1'^^]^^ • For example, ici{q±;q±) = -iici{q±] q±) ■ 
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Of course, {i^iC)ci{(l±] (1±) = ^ci{(l±', (1±) ■ (^'CC)c/(?±; fe); because the two operators of our 
representation commute each other. Then, (ci{(l±', (1±) = 2(i^C)c«(9'±; Q'±) ■ ^ci{(l±](l±) and 
Cc/(q'±; (1±) = ~Kci{.<l±] (1±)- The dependence on time of these random variables can now be 
obtained from the equations of motion (jHHjl . (jTTj) . 

Once we have stated the rules to define the random variables on configuration space 
associated to each operator we can pick E,ci and as the two independent variables, noticing 
then that {i^()d = \Cd/ici- In the vacuum the four equally probable random events in the 
configuration space correspond to the four possibilities: = Ccz = 

The whole formalism can be repeated for any quantum state \^) other than the vacuum. 
First, we need to expand the state |\E') in the basis |\1') = ''p{+]+)\q+'A+) + 

^(+; ~) k+S 9-) + V"!"! 0'+) + "^("S Each PGP correspond to one of the 

points g^) in the configuration space and its probability to happen randomly is equal to 
the modulus squared of the amplitude of the wavefunction, |(g|\E')p. The four events are no 
longer necessarily equally probable, but the classical values on each of them of the random 
variables assigned to the operators that form our complete representation of commuting 
observables are still given by their corresponding eigenvalues fl43p - (Pmi . Then, the action of 
the operator 0{t) on the state |\E') is written as the action on the vacuum of some operator 

g{t-l,^MO-- 0{tm = (xo;o(t)-l + Xl;o(t)-^ + Xl;lW-(^^C)+XO;l(t)-e(^eC))|0;0) ^ 

Q(t; l,^,i^()\0; 0). This step is always possible and uniquely defined because the action on 
the vacuum of each one of the four operators and ^(i^C)) excites a different mode in 

the basis of eigenstates of the hamiltonian (f^H|) . Finally, the relation between the states 
|\E') and |0;0) is inverted |0;0) = J'{l,^,iC,C)\'^) and introduced in the previous expression 
to get the desired result: 0(t)|^) = g{t; 1, if'C) ^^(1, ^, ^I^C)I^) = Vfit; 1, if'OI^)- This 
last expression can be used to define the classical value of the operator O in |g) when the 
quantum system is described in its state |\1') as Odit; q) = Vf{t; l,^ci{q), {i(,Od{q))- 

We see that although each PGP corresponds to one of the eigenstates the way 

how we actually define the classical values of physical observables explicitly depends on 
the quantum state of the system |\1'). Therefore, the PGP cannot be identified with the 
quantum state described by the eigenstate Instead the whole set of four random PGP's, 
which according to ((221)- (EH) should still obey classical equation of motion, does describe the 
dynamics of the system in its quantum state |\E'). The state of the system, therefore, enters 
this description by: first, fixing the different probabilities of each of the four points of the 
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configuration space; second, defining the initial conditions on the set of classical equations 
finn|) - (jMj) for each of the PCP's; third, defining the algebraic relations between the values of 
different physical observables on each PCP. 

The formalism can now be straightforward generalized to a fermionic system containing 
many modes. We define a couple of hermitic operators C,i,^i associated to each one of these 
modes. The integer index / = 0, 1, 2, ... labels the different modes. The set of operators 
i^o^i, i^i^2,--- forms a complete representation of commuting observables and the Hilbert 
space of the states of the system can be linearly spanned in the basis of common eigenstates 
to all the operators in this representation. Each common eigenstate can be labeled by a 
sequence of + or — signs (for example, |+,+,— ,+,—,— ,—,...)) each input corresponding 
to one of the fermionic modes. In total, 2^ different eigenstates, where Af is the number of 
different modes. A + sign in the zeroth position corresponds to an eigenstate with eigenvalue 
^ for the operator ^o- Inversely, a — sign correspond to an state with eigenvalue ^ for this 
operator. And similarly, a ± sign in the l-th position corresponds to an state with eigenvalue 
^ for the operator I > 1. In the quantum state |\t') = '4^{(1±^)\(1±^) each of the 

2^ random pseudoclassical paths |g) which describe the system has a probability to 
happen. In the vacuum all the paths are equally probable and, therefore, this probability is 



Let us now turn back to equations ()23|) and ()24p . which describe the statistical behaviour 
of the random variables over the ensemble of PCP's. They relate the first and second 
momentum of the random variable Od to the expected values of the operator O and O'^O, 
respectively. The proof of these equations that we gave in the previous section also proves the 
assertions for systems with many fermionic modes. We noticed then that ()24|) holds only on 
average over the whole ensemble of PCP's but does not necessarily hold on each specific PCP 
and understood that, as a consequence of the non-commutative relations between quantum 
operators, the usual algebraic relations between classical observables are not necessarily 
fulfilled by the time-dependent values Oci{t) defined on the non- interfering paths. There are, 
nevertheless, a certain class of operators, which we call "collective" operators for reasons 
that will become clear in the next section, for which (O^O)c/ — {Oci)*Oci ~ not only in the 
average sense of equation (|^ . but under the stronger requirement that 



1/2^. 





(48) 
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In this sense (0'''0)ci ~ iOci)*Oci on each PCP independently of the quantum state that 
they describe and we can conclude that values on paths of collective observables approxi- 
mately recover the usual algebraic relations between classical observables. This condition 
is naturally satisfied for some operators which depend on the dynamics of many modes 
because for these operators the random variable (0''"0)c/ — {Oci)*Oci can be expressed as 
a quadratic function of many other independent random variables associated to individual 
modes and, according to the central limit theorem, condition ()48|) is satisfied. Moreover, 
higher statistical moments of the collective random variable are highly suppressed. 



IV. PCP'S IN FERMIONIC FIELD THEORIES 



We are now ready to develop the formalism of PCP's for the Dirac theory of free spinors. 
We consider a single fermionic field ip{x) in Minkowski space-time, whose lagrangian density 
is Q 

C = ij (id^Y - M) ijj, (49) 

where j'^ are 4D Clifford matrices and M is the mass of the spinor field. The hamiltonian 
of the system is: 

n = j (fx V^(x) {-i^^dj + m) ^(f). (50) 

In order to regularize the theory in the infrared limit we impose periodic boundary con- 
ditions on the large three-dimensional box [0, X] x [0, X] x [0, X]. The spinor field ip{x) and 
its conjugate il){x) = ip^x)'')^ can then be expanded in a discrete series of Fourier modes: 

+N 1 

m= E ^ E + ^^k^e"'™^'') ' (51) 

n\,n2,nz=—N s 
+N 1 

m = E ^ E {ajnte-'-'''-'^'' + htvte^'-^--!^) . (52) 

ni,n2,nz=—N s 

The index s takes two possible values which correspond to two possible helicity ori- 
entations for the fermionic mode labeled by the index n = (ni, ^2, jt-s). The infinite 
series of Fourier modes has been cut by the ultraviolet regulator N/X. The opera- 
tors a^, and their hermitic conjugates al^, 61^ obey canonical anticommutation rela- 
tion: {a^,a^^} = {bfi,b^} = d^'^'darn and any other anticommutator between these op- 
erators is equal to zero. They create and annihilate quanta of the fermion and its an- 
tifermion, respectively. The Dirac spinors are two linearly independent solutions to 
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the equation [—i'yWj + M]u^e~^'^'^^'^'^^^ = K^'y^u^e'^'^^^^'^^^ , while v'^ solve the equation 
[-i7^'o»j + M] t;^e-2™-^/^ = -fi;^7%ee-2™-^^A, where = [|^p + M2]^^^ They obey 
orthogonality relations of the kind u^Jv^^^ = f ^mI'^ = and are normalized such that 

Once we put ()5ip and ()52|) into (jSUj) we obtain the following expression for the hamiltonian 
of the system 

^ = E '^n^ - ^) + iK% - ^)) , (53) 

n.s 

which describes a collection of free fermionic modes. 

We introduce for each mode, labeled by the pair of indexes s, n, a couple of hermitic 
operators and as we did in the previous section. The task of choosing a complete 
representation of commuting observables is much simplified after having regularized the 
theory in the IR and UV: then we have a large but finite number of fermionic modes that 
we can order in a sequence. The set of hermitic operators '^^o^i5'^^i^2, ...forms a complete 
representation of commuting observables and we can expand any quantum state of the 
system in the basis of common eigenstates \q) to the operators in this representation. 

Following the formalism that we have developed in the two previous sections we can give 
to each physical observable 0{x,t) a "classical" value Oci{x,t) at each of these common 
eigenstates which are so promoted to describe space-time dependent pseudoclassical 
field configurations (PCFC's) whose probabilities to randomly happen are, again, given by 
|(g|\l')p. We remark again that although each PCFC corresponds to a certain eigenstate 
\q), the way how we actually assign classical values to the physical observables, exploiting 
identities like (fTTj) . explicitly depends on the quantum state of the system and, therefore, 
the path cannot be identified with the quantum state described by |g). 

The state |\1') of the quantum field is then represented as a linear combination of randomly 
distributed PCFC's. Each field operator, in particular ip^x) and its conjugate ijj{x), is 
assigned a time-dependent value on each of these paths. The operator equation 

[td^r-M]^{t,x) = 0, (54) 



and its hermitic conjugate 

^(t,f)[^9^7'^ + M] =0, 



(55) 
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when are applied on the quantum state of the system 

[id^Y - M]ilj{t,x)\^) = 0, [z9^7^ + M] 1^) = 0, (56) 

guarantee that on each of the random PCFC's in the set that describes the quantum state 
the classical values ip{t,x)ci and ip(t,x)ci obey classical equations of motion: 

[td^^>' - M] ^(t, x),i = 0, ^P{t, x)ci [id^Y + M] = 0. (57) 

On the other hand, the initial conditions on each of the PCFC in the set that describes 
the state, as well as their probabilities to randomly happen, are different for each quantum 
state of the system, as we have discussed in the previous sections. Let us remark that 
{4')ci = (V'W°)c/ is determined from the action of the operator on the ket-state ipl"^) and, 
therefore, is not equal to {ipdYl^ which results from its action on the bra-state 

We can now consider the addition of a linear interaction term to the quantum field free 
theory that we are discussing. If the coupling is weak enough we can work in perturbation 
theory. Then, the Hilbert space of states can still be built as the Fock space of free quanta and 
the algebra of the quantum operators is not altered. So, the whole formalism of PCFC's can 
be repeated. Now there is a new interacting term that appears in the hamiltonian of the sys- 
tem and modify the dynamics of the quantum operators, 0{t) = Q^(iio+^Hi„t)tQ^-'i.iHo+XH,„t)t^ 
in the Heisenberg picture. In consequence, it also modifies the time- dependent values of the 
field operators on each PCFC. We can see from the operator equations ()54|1 and that 
ip{t,x)ci and ip{t,x)ci will still obey exactly classical equations of motion because the term 
added to the equations is linear in the operators ip(t,x) and ip(t,x). 

As a simple example, we can add a Yukawa coupling term to the lagrangian ()49|) . This 
interaction adds a new term to the operator equation: 

[id^Y - M] ^{t, x) + g(f){t, x)^{t, x) = 0. (58) 

This equation implies, when it is applied on the quantum state that the dynamics of 
the pseudoclassical values on its PCFC's gets modified: 

[td^Y - M] ^/j{t, x),i + g(j){t, x)cMt, f)e/ = (59) 

and a similar equation can be obtained for {ip)ci- In the new term we have made use of 
the identity {4'ip)ci = {(p)d{4')ch which can be easily proved because in a weakly interacting 
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theory the operators and ijj act on different sectors of the Hilbert space, and have taken 
advantage of the hnearity of (j^Hj) in the fermion field operator. The interaction term does 
not modify the initial conditions to be fixed on the PCFC's that describe the quantum state, 
which depend only on the relations between operators when they act on the quantum state. 
Neither it does modify the random probability of each path. A discussion of, in general, 
non-linear interactions will be presented in a forthcoming work. 

In the picture that we have presented quantum states of fermionic fields are represented 
as linear combination of randomly distributed paths that obey classical Dirac equations and 
do not interfere between themselves. The formalism is reminiscent of the consistent histories 
approach of Gell-Mann and Hartle Q], wherein the Feynman paths [l^ which contribute 
to a certain process are grouped into equivalence classes or non-detailed histories that are 
approximately incoherent. On the contrary, in the formalism that we have developed the 
random paths do not interfere at all and obey "deterministic" equations of motion at the 
price of redefining the algebraic relations between the values on paths of physical observables. 

We can use this formalism of PCFC's to describe quantum fluctuations of fermionic fields 
in terms of pseudoclassical stochastic proccesses. In particular, we can use it to describe 
local fluctuations of globally conserved numbers using classical concepts. We present here 
two examples that could be interesting: local fluctuations of the energy and fermion number 
in a finite volume V contained in the whole box [0,X]^ and parametrically larger than the 
UV cutoff that we have introduced. 



x) which describes the energy contained 
13 and discussed in M in the context 



The operator Hv = Jv i'ix) {—i'-yWj + M) ifj 
in the finite volume V was already introduced in 
of a bosonic field theory. Similarly the fermion number contained in the same volume is 
described by the operator By = fyd^^ ^/'(x)7°?/'(x). 

If we want to find the "classical" expressions for these two operators in the set of PCFC's 
that describe a generic quantum field theoretic state we would need to write the action 
on this state of each of the operators in terms of the action on the same state of some function 
of the operators in a complete representation of communting observables. We have done this 
excersise in previous sections and, therefore, we will present here a different approach that 
will shortcut the way towards this aim and will give us a different perspective on PCFC's. 
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We know from equation ()24|) that 

< {nv{t))cl >PCFC = f d^X < [-i-f^dj + M) ^{t,x))cl >PCFC 

= j d^'x < (V'd(^,^))^7° (-11' dj + m) ^d{t,x) >PCFC, 

and similarly 

< iBvit))^l>PCFC = [ < {lp{t,x)l°'^{t,x))d>PCFC 

Jv 

= / d^X < {lpd{t,x)yipd{t,x) >PCFC ■ 

Jv 

Although these equations hold only on average over the whole ensemble of PCFC's and not 
necessarily on each of them, we gave in ()48|1 a condition that ensures that 

{nv{t))ci ^ f d^x (Mt, {-^l'^, + m) ^,i{t, x), (60) 

and 

(ey(t)),,~ / rf^X {Mt,x)yiJd{t,x), (61) 



on each PCFC, if the random variables 
and 



{nv)ci - Iv d'x {Mx))h' Hl'd, + M) i;,i{x)_ 
{Bv)ci — Jv {4'd{x)yipci{x) are quadratic in many independent random variables. 
We called this kind of observables collective observables and will show below that these two 
operators are in fact " collective" operators, but for now let us just assume it to get to some 
conclussions: 

a) We can describe the dynamics on paths of " collective" observables like (jUUI) and (jUTj) 
in terms of classical stochastic concepts, from the classical Dirac equations (j37|) or (fS^ 
that ipd{t,x) obeys. Only the initial conditions and the actual distribution of paths are 
constrained by the quantum state of the system |\E'). 

b) Then, the time correlation function of any operator can be defined on each PCFC as 
usual in statistichal mechanics [l^, fo{ti — ^2) = i0d{ti))*0d{t2) — {Od{t2))* Od{ti) , and 
the lifetime of the fluctuation is defined as the inverse width of its Fourier transform. We 
proved in that the average over the ensemble of PCFC's coincides, as it should do, with 
the usual definition of the time scale of quantum fluctuations. 

We will conclude this section giving a formal proof that justifies fl6()|l and ()61|) . Let start 
introducing the Fourier expansions (j3T|) and ()52|) in the expressions for the operators Tiy 
and By- Their classical values on paths can then be expressed as [fsj ]: 

inv)d = E E + &-n^-n)(47°4 - &-!^7%-'^)l^^ Fin -m,X), (62) 



n,m s,s' 
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ra,m s,s' 



On the other hand the approximate expressions (jUm) . ()61|) can be expressed using the same 
Fourier expansions as: 



{ny):r'^ = EE'^-i [&i]:.^^^)([4]c/7°4 - kIIc^tV^)] F{n-m,x), 

n,m s,s' 

(64) 

{Byrr"^ = E E [([^TcfiH + [b%]:iV-'.MaiWui + rlWv!:^)] F{n -m,X). (65) 

n,rn s,s' 

We need to prove that {Hv)ci — {'Hv)'d^^°^ and {Bv)ci — {Bv)ci^^°^ are bihnear in the random 
variables associated to the single modes. Comparing both expressions we can see that it is 
enough to prove that expressions of the kind [a^Ja^]ci — [0'fi]ci[(^^]ci, which involve only two 
different fermion modes are bilinear. We know from that they are. 



V. SUMMARY 



We have extended to systems of fermions the formalism of pseudoclassical paths that we 
recently developed for systems of weakly interacting bosons and have shown that fermionic 
quantum states can also be represented in the Heisenberg picture as linear combinations 
of randomly distributed paths which do not interfere between themselves. Every physical 
observable was assigned a time-dependent value on each path in a way that respects the 
anticommutation relations between fermionic operators and, in consequence, these values on 
paths do not neccesarily respect the usual algebraic relations between classical observables. 
Each path depicts a collection of pseudoclassical harmonic oscillators with constrained initial 
conditions. 

We used these paths to define the dynamics of quantum fluctuations in systems of fermions 
without reference to an environment or any additional external system. Then, we selected 
"collective" observables which depend on many fermionic degrees of freedom and found that, 
as we found in for collective operators which depend on many bosonic modes, their values 
on PCP's do approximately fulfill the usual algebraic relations between classical observables 
and, therefore, the dynamics of these collective fluctuations can be described in terms of 
unconstrained classical stochastic processes. 

In this setup, we showed that quantum fluctuations of fermionic fields obey classical Dirac 
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equations and described the dynamics of local fluctuations of globally conserved fermion 
numbers. 
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